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Abstract 

We propose a new interpretation of the BFSS large N matrix quantum mechanics anal- 
ogous to a novel interpretation of the IKKT matrix model where infinitely large N matrices 
act as differential operators in a curved space. In this picture, the Schrodinger equation in the 
BFSS model is regarded as the Wheeler-DeWitt equation which determines the wave func- 
tion of universe. An explicit solution of wave function is studied in a simple two-dimensional 
minisuperspace model. 
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1 Introduction 

In matrix model proposals for non-perturbative definition of M-theory (the BFSS model flj) and 
superstring theory (the IKKT model [2J), realization of quantum theory of gravity has always 
been a major concern. String theory in its perturbative formulation contains graviton excitation 
and it gives a consistent quantum theory of gravity. M-theory is defined as the strong coupling 
limit of type IIA string theory and should include the eleven dimensional supergravity as low 
energy effective theory. It has been expected that their quantum aspects can be captured in 
terms of the IKKT and the BFSS matrix models. 

However, it is not easy to identify these matrix models with quantum gravity, because 
general covariance is not manifest in their formulations. Rather, they are formulated as large 
N limit of gauge theories; namely dimensional reductions of AA = 1 super Yang-Mills theory in 
ten dimensions down to 0-|-0 dimension (IKKT) or 0-|-l dimension (BFSS). Traditionally, these 
theories have been interpreted as a systems of infinitely many D-branes. In such interpretation, 
the gravitational effect appears rather indirectly through the loop effect of open strings stretched 
between each D-branes. Such idea has already been examined, and it was found that the loop 
effect reproduces the graviton exchange between two D-objects [1112]. The graviton exchange 
was extensively investigated further in higher orders of loop expansion and derivative terms, and 
correspondence between the matrix theory and supergravity was reported OH]. As another 
attempt the induced action of gravity by a loop effect in the matrix model is discussed in [5] . It 
has been a difficult task, however, to find a way of realizing quantum gravity in terms of matrix 
model without relying on any perturbatioiu. 

In recent years, an alternative identification of gravitational degrees of freedom in matrix 
model have been proposed by Hanada, Kawai and Kimura [8]. These authors showed that a 
matrix can be identified with a covariant derivative in a curved spacetime in large A^ limit. Once 
we apply this identification to the IKKT model, for instance, we find that the matrix equation 
of motion gives Einstein equation in vacuum. In this paper we will call such an identification 
in matrix model as "HKK interpretation" for convenience. Following this approach, we do not 
need any loop calculation, and the Einstein equation is derived quite straightforwardly. Such 
feature is quite attractive even though there is still subtlety of regularization due to infinite size 
of matrices. 

Past studies along this line have been devoted to the IKKT matrix model [9l \TU[ \TT\ . In this 
paper we address an application of the HKK interpretation to the BFSS matrix modelcl This 
academic application immediately gives us some interesting results. We recall that the BFSS 
matrix model is a quantum mechanics whose dynamical variable are time-dependent matrices. 
It is, therefore, natural to expect that the Schrodinger equation describes quantum evolution 
of spacetime. This is remarkable advantage of our formalism in contrast to the IKKT matrix 
model whose matrices are time independent. The Schrodinger equation becomes in the HKK 



^ As for recent studies on string theory aspects of the models, a hght-cone superstring action has been derived 
from IIB matrix model 16 , and vertex operators and scattering amplitudes of superstring are reproduced in ^. 

^Although we study the large A*' (bosonic) Yang-Mills quantum mechanics with d matrices, still we call it the 
"BFSS" model. Similarly we call the large A'' Yang-Mills reduced model with d+ 1 matrices the "IKKT" model. 



prescription: 

i-^it, X) = H{x, d/dx)-^{t, X) -^ i-^{t, V) = H(y, d/dv)^{t, v), 

where V is a covariant derivative in a curved space described by the vielbein and the spin 
connection. A wave function obtained by solving the above equation describes quantum evolution 
of the spacetime. We will also address the meaning of expressions such as ^ = ^{t, V) and d/dV 
which look somewhat symbolic. 

This paper is organized as follows. In section [2l after reviewing the HKK interpretation 
and the quantization procedure of the BFSS model applicable to finite A^ case, we apply the 
interpretation to the BFSS model and obtain a Hamiltonian in terms of the HKK interpretation. 
The Schrodinger equation is derived straightforwardly from the Hamiltonian after quantization. 
This process is in contrast to traditional approaches of quantum gravity where the Schrodinger 
equation (the Wheeler-DeWitt equation) is obtained as a Hamiltonian constraint. In section [3l 
we examine a minisuperspace model in two dimensions as an explicit example. We find that 
the time independent Schrodinger equation can be solved analytically. We also obtain a time 
dependent wave function by constructing a wave packet. Expectation values and dispersions 
of various geometrical quantities such as metric and torsion are evaluated. Finally we discuss 
physical implication of our result. 

2 HKK interpretation of BFSS model 

2.1 HKK interpretation for large A^ matrices 

We begin with a brief review on the HKK interpretation of matrix models proposed by [8j 
where the authors argued how to interpret large A^ matrices as covariant derivatives in curved 
spacetime. Their formalism can be applied to arbitrary mainfold at least formally, and a case 
of some compact manifolds had been studied in detail[TO]. In their proposal, a large A^ matrix 
is identified with an operator which acts on a section of vector bundle Ereg over D{= 1 + d) 
dimensional curved manifold M endowed with structure group G, where G = Spin{l,d) is local 
Lorentz groupO. The relation between large A^ matrix X/a) ^^'^ the covariant derivative V^ on 
a manifold M is given as 

X(„) = iR^^fig-')Va. (2.1) 

The matrix R,^'^{g~^) is the vector representation of local Lorentz group G {g ^ G). While the 
index a transforms as Lorentz vector with respect to the action of G on T{Ereg), the other index 
(o) remains unchanged which is a index of global SO{D). This feature enables us to identify 
(j2.ip as an endomorphism on T{Ereg)^ i-e., matrix acts on this space. In this formalism, the 



^ The section of a fiber bundle V{Ereg) is a set of smooth maps from a coordinate patch of M to some vector 
space Vreg- Vreg ~ {f '■ G ^ C} is reducible and has following decomposition 

Vreg ^ 0r ( v r ' ' ' Vr ) : 

" V ' 

d,. 

where Vr is the space of irreducible representation r of G, and dr is its dimension. 



covariant derivative V^ is interpreted as a linear map from T(Ereg) to TM T{Ereg), where 
TM is the tangent bundle on M. It can be written explicitly as 



Va = e,^(a^ + V^a,), (2.2) 

where e and to are vielbein and spin-connection respectively, and O is Lorentz generator whose 
explicit form depends on a representation on which it acts. 

2.2 BFSS - matrix quantum mechanics 

Before applying the HKK prescription to the BFSS model, let us recall quantization of the 
model in a way that is convenient for later argument. The bosonic part of the BFSS action in 
Lorentzian signature is given by 

Sbfss = j dtL = Jdttr (^-^[Do,X^,)][D^,X^^^] + i[X(,),X(,)][xW, X(^-)]) , (2.3) 

where Dq = dt + iAq is the gauge covariant derivative. The indices (i), {j) run only in d spatial 
directions. We shall work in the ^o = gauge such that [L'o,-'^] = dtX. We employ the matrix 
notation as 

X^l(t)=x^2(t)t^,, (2.4) 

where t^^ is the generator of SU{N) algebra in adjoint representation so that a and b run from 
1 to A^ and indices A from 1 to N"^ — 1. x^^ is the degrees of freedom and t^^ are the basis. The 
conjugate momentum of Xu\ is defined as 

(i) _ dL _ • (i) 

(jX[i)ab 

where X denotes dtX. The Legendre transformation gives us the BFSS Hamiltonian 

H = tr|in(,)n«-i[X(,),X(,.)][xW,x(^-)]|. (2.6) 

The gauge constraint arising from our choice of Aq = gauge is 

[X»,n(,)]=0. (2.7) 

The Poisson bracket of Xu\ and IT'--'' is written as 

^Wa.'n2}^^ = ^^i)^'^hABtiC, (2.8) 

where h^B is the inverse of the metric h = t^^t^. Upon quantization the Poisson bracket 
becomes the canonical commutator, and n^*-* acts as —id/dXij\ on the wave function '^{X). 
The Schrodinger equation then becomes 

Having reviewed the matrix quantum mechanics, we will apply the HKK prescription to the 
above argument in the next subsection. 



2.3 Applying HKK to BFSS Shrodinger equation 

We set the matrix Xu\ (t) as a covariant derivative: 

X(,) = iR^,)'Vi = iR^ifie/di + LO^'^Ojk), uji^'' = e/iOi^\ (2.10) 

where Ru\^ {g"^) belongs to G = Spin{d) vector representation and all the indices run only in 
d spatial dimensions. Note that the expression (j2.1Up is analogous to (|2.4p if we regard dj and 
Ojk as basis, and the vielbein and the spin connection as degrees of freedom. We shall make 
more comments on this formal argument in Appendix |Al We also require each index in (j2.10p 
to transform appropriately in d dimensions. Thus the covariant derivative in (j2.10p is associated 
with a d dimensional manifold for each t. One may wonder whether this manifold corresponds 
to a time slice of some d + 1 dimensional manifold whose timelike Killing vector is given by dt ■ 
Indeed, in Appendix|Bl we demonstrate that one can choose a local frame equipped with metric 
—dt^ + h^{t,x)ijdxdx such that Xu\{t) and Dq are embedded into the d+\ dimensional IKKT 
model, although we do not pursue along this line of reasoning here. 
The commutator becomes|j 



with 



where 



[X(i),X(,)] = -i2(,)*i2(,/-[V„V,], (2.11) 



[V^,V,] = [e/V7,e/Vj]=ri/ai^ + i?,fOfcz, (2.12) 



K — ^La „K , , . K 



e[f5ie/ +.;[,/, (2.13) 

p Kt — fy kl , rp K , . k 
JXij = l-Cij -\- lij LOk 

Here T^ is the torsion and TZij ' is the Riemann curvature tensor. The reader may refer to 
Appendix [O for derivation. If the torsion free condition Tij = is satisfied, then a and cui^ 
are no longer independent of each other. We have learnt that the first proposal discussed in 
[8] considered torsion free case for simplicity. In general, however, matrix configurations do not 
necessary to satisfy the torsion free condition. We remark that torsion in the IKKT model has 
been discussed in pT] . 

Next we apply the HKK prescription to the conjugate momentum (12. Sp . We find 



n» = xii) = i6^')(^)R^i^\ei^di + u;.^''Ojk)- (2.15) 

In ()2.10p and ()2.15p . only e/ and LOi^ and their time derivatives are dynamical, while dj and 
Ojk are independent of the choice of geometry. Therefore it is convenient to introduce a phase 
space defined by {ej^, coi^'^} and their conjugate momenta 

K/ = e/, ^-f =c:;,^-^}, (2.16) 



''The representation matrix R(i)'' {g ^) can be moved into left side of tfie covariant derivative Vi. It was 
explained in [8] in detail. 



where the Poisson brackets are defined by 

{e/(x),7rY(2/)}pB = d,,6''6''{x - y), {u;P\x),7r'^riy)}pB = SuS^'^S^^d^x - y). (2.17) 
Then using the above definition, the Poisson bracket between Xu\ and IIq) can be evaluated as 

[x^^{x)^pM'\yU} ^^ = %)^'^'J''(^ - y) [-did'i^pi^s - {Oim)o.p{o"^)-,s)) , (2.18) 

where we write matrix indices along fiber direction explicitly in Greek subscripts. We have also 



used the orthogonal property Ru) Rtj) = S(i)(j)- Note that there is a "factor" in (j2.18p which 



multiplies SlJd {x — y). We also learnt that similar factor also appears in ()2.8p . In fact, such 
analogy between (|2.18|) and (|2.8|) can be seen by identifying di and Oij as t^^ in (|2.4|) . but we 
can not put forward this analogy further since t is a matrix but di and Oij are not. Parts of 
these combining with Ru\'' may provide corresponding basis. See Appendix [Al 

Having defined the phase space, we are ready to construct Hamiltonian which governs dy- 
namics of our model. It can be obtained by applying (|2.1U|) and (|2.15|) to the BFSS Hamiltonian 
defined by (j2.6p . The trace in (j2.6p should be performed over a complete set of function in 
C'^{Eprin)- In this way we have 

H = - j d'^x^,tr^g^-{^t^dK)K^dL) + \{7^r'0,k){^r''^0im) 

+i(r,/a,,)(ri/9L) + -^{Ri,'w,,){R,,'''Oki)\ , (2.19) 

where we have used (JB.IOP for tr in the large A^ limit. The operation tvg is defined as trgF = 
J dg{x,g\F\x,g) with a Haar measure dg for g. The matrix -R{i)* has disappeared thanks to its 
orthogonality. There is no cross term because trgO = 0. There are the trace operations such as 
J d^xd^ and tr^O"^ and they formally diverge due to tracing over a regular representation with 
infinite dimensions. Here we simply assume that we have employed a suitable regularization 
procedure, which could be the heat kernel regularization or elsewhere suggested in |10j . 

Having obtained classical Hamiltonian, we can now quantize the system by replacing the 
Poisson brackets in (j2.17p with commutators. We shall work in a representation such that e and 
(jj becomes diagonal. With this choice, vr^^ ^ and vr^ -^ are promoted into operators, 

<' = -^/7, <'' = -^-^,. (2.20) 

Note that there is a problem of operator ordering due to the presence of the determinant factor 
in the Hamiltonian. Choosing an ordering prescription such that the Hamiltonian preserves 
hermiticity, the Schrodinger equation (j2.9p becomes 



4^ = /^'^v^*^4^^^^^-)^(4^-) + ^(^^^•^)(4^^-) 

;(T,/a^)v^(r,/9i) - liRij''OM){Rij'^''0^n) W. (2.21) 



1 1 



In addition, the gauge constraint (|2.7p . while imposing on the wave function, takes the fohowing 
form 



{T^dK + R''Oij)^ = 0, 



(2.22) 



where 



nK 



R'^ 



\dLet) 



\dLe. 



K' 



K 






e,-dK^l^-e,-dKu:l=+urd:, 



■ I K 

ik, kj 



UJi UJ^' + iVii UJ^ 



k,\ «i 



^n ^k 



(2.23) 



The equation (j2.2ip together with the constraint (j2.22p dictate the evolution of quantum system 
with degrees of freedom given by vielbein e and spin connection to, which in turn determine 
gravity in the classical level. To summarize, we regard that the Schrodinger equation (j2.2ip and 
the gauge constraint (|2.22p are the equations of quantum gravity which is realized in the BFSS 
quantum mechanics at large N limit. 

In the case the Hamiltonian in ()2.2ip is time independent, by substituting ^(t, x) = e~ '^ e{x) 
we obtain the time independent Schrodinger equation 



d x^Jhd tr 



1 1 



^2'6ei^ 



dK)V^i^-r^L) + 



6ei^ ' l^bbJiJk 



'ak ' 



-O 



Im) 



1 1 



/id 4 



-Mj'^dKWhdiTij'^dL 






-{Rij Okl){Rij 



^E = E-^E- (2.24) 



The appearance of this equation, at first sight, is analogous to the Wheeler-DeWitt equation in 
canonical quantum gravity [12j. One might, however, question their similarity mainly because 
the Wheeler-DeWitt equation is seen as a Hamiltonian constraint without time evolution due 
to reparametrization invariance. Nevertheless it is known that the Schrodinger equation can be 
rewritten as a Hamiltonian constraint through time reparametrization, i.e. t = t{T) [I3j. In this 
sense the equation (j2.2ip and (j2.24p indeed play the same role as the Wheeler-DeWitt equation 
does. Here we briefly recall the argument. We introduce the lapse function N{t) such that 



dt{T) = N{T)dT 



(2.25) 



in the action (j2.3p . The variable ^(t) is regarded as an independent variable with equal footing 
as X{t). Introducing a Lagrange multiplier vr which plays the role of the conjugate momentum 
with respect to t(r), we have 



S = Tr dr 



1 



2N{t) 



dr 



dt{T) 



dr 



N{r] 



(2.26) 



The conjugate momentum with respect to X^'^\t) becomes 

1 dX^^ 



p(i) 



N dr 



then we obtain 



S = Tr dr 



P(,)XW+^t-A^C, 



H 



(2.27) 



(2.28) 



where 

CH = \pfi) + \[X^\x'^^)f + 7.. (2.29) 

In this way we obtain the canonical system with variables (X, t) whose conjugate momenta are 
(P, tt) respectively and the Hamiltonian is given by NCh- We obtain a constraint Ch = by 
the variation of A^. Clearly the action is invariant under the reparametrization t ^>- t' = t'{t) 
in which A^ -^ N'(t') = dr/dr'N^r), and all the others transform as scalars. The Lagrange 
multiplier vr plays a role of E in (j2.24p . 

Now one can consider the WKB approximation as a semi-classical limit of the large N BFSS 
model. In the WKB approximation, a wave function is governed by a classical action evaluated 
at saddle points, which satisfy the classical equations of motion. In the present case we have 

[Z)°,[Do,^«]] + [^(^'\[^0-),^«]] = 0, (2.30) 

[X«,[X(,),Do]] = 0. (2.31) 

By using the HKK interpretation, they become (see Appendix [D] for derivation) 

e^' + {V'n^)e/ + 7J"'Tk^"'TlJ-T/u;)J + Ri' = 0, (2.32) 
^iki - Ci ^e/u;jM + T^mi"{RZki " r/^^rp^^^^.H) + V^{Rjiki - Ti.'^uj^ki) = 0. (2.33) 

Hence the semi-classical limit of the large A^ BFSS model describes a theory of gravity with/without 
torsion, even though the Hamiltonian itself looks very different from that in the Einstein the- 
ory. In other words, quantum nature of the BFSS model could be different from what has been 
derived from quantization of the ordinary Einstein-Hilbert action. However, we should remind 
that above consideration is somewhat naive therefore has to be refined at least with respect 
to the following points: the first point is that although the form of equation keeps the form 
of Einstein equation, general covariance is lost in our application of HKK to the BFSS model, 
namely complete classical Einstein equation appears only in spatial directions. The second point 
is that in the quantum mechanics the time evolution of expectation value of an operator A obeys 
Hamilton's equation, 

j^{A)=i{[H,A]). (2.34) 

Then expectation values of the position operator x and the momentum operator p in quantum 
mechanics obey classical equation of motion. Here we have additional metric determinant in 
the large N Hamiltonian and its time evolution makes time dependence of expectation values 
depart from classical matrix equations of motion. These discrepancies from classical Einstein 
theory will be significant where spacetime is far from the flat one without torsion. 

3 Minisuperspace model 

In this section, we will examine our new interpretation in a minisuperspace model as an explicit 
example. To be precise, we will consider an exactly calculable toy model of two dimensional 



universe to clarify our proposal given in the previous section. We start with a simple (gauged) 
one matrix model without potential term: 



£ = hr[Do,Xf. 



(3.1) 



It seems to be almost free, but the metric determinant factor which appears after taking large 
A*" limit gives nontrivial dynamics. Assuming our toy model of universe is given by the following 
metric: 



ds^ 



-dt^ + 



1 



yity 



:dx' 



(3.2) 



This is a two dimensional version of the Robertson- Walker universe. In the assumption here the 
function y{t) only depends on time. It makes our setting to be a tractable problem. Following 
the HKK interpretation, we obtain the covariant derivative: 



Do = ^§-^, X = ^yit)^, 



(3.3) 



as well as the time independent Schrodinger equation: 



A — - — 

dyydy 



A 



tvg I dx{dxf. 



(3.4) 



Here A is a positive and divergent quantity which is expected to be regularized by some reg- 
ularization procedure like the heat kernel regularization proposed in p^, and in the following 
we will treat A as if it is a finite quantity. The factor 1/y sitting between the derivatives is 
originated from the metric determinant factor in the Hamiltonian. The gauge constraint ([27 
on the other hand is trivially satisfied in this setting. From (|3.4p . we obtain 



d 1 d 
dyydy 



+ e 



^ = 0, e = E/A. 



(3.5) 



A solution to (13.51) for e > n is found to be 



i: 



*. 



yJ+1 ( T7\/ey2 



V3 



(3.6) 



where Ju{z) is the Bessel function. The wave function ^^ is orthogonal 

^e{y)^e'{y)dy = 5{e-e'). 



(3.7) 



We have provided a proof for the orthogonality in appendix [El 

Using this orthogonal basis (j3.6p . time dependent wave packet is constructed as follows: 



^(i,y) 



de C(e)e-*^^'^,(y), 



(3.8) 







^It corresponds to plane wave type solution in the BFSS quantum mechanics. 



where C(e) is arbitrary function that is normaUzable. Here we take the following Boltzmann 
form as an example 



C{e) 



r(5/3) 



€3e 



-/3e 



(3. 



such that normalization condition f deC{e)'^ = 1 is satisfied. /? is a parameter to characterize 
the wave packet. Then equation ()3.8p becomes 



^{t,y) 



(2/?)i 



y 



r(5/3)3i(^ + ,^i): 



■ exp 



y 



9{l3 + iAt) 



(3.10) 



The wave packet becomes zero as t — > itcxD or y ^ oo and its peak is located at {t, y) = 
(0,(15/3)^''^). We observe that the wave packet exponentially decays for large y although the 
metric ()3.2p shows a singularity at y — > oo. However there is no pathology about this because 
the wave function itself is regular and smooth for all values of y. 

The time evolution can be seen by evaluating the following expectation value: 



(y"") 



dy^*{t,y)y"'^it,y) 



r(f) 



2+2 



/3" + A H 



(3.11) 



where we have used (|3.10p . One can easily calculate various expectation value of operators based 
on (j3.1ip . To our most concern, the scale factor is found to be 



( 9xx ) = {y ) 



m 



2 
2\ 3 



P 



/32 + AH"^ 



(3.12) 



We observe that around t = the expectation value of size of universe keeps a constant value, 
and late time behavior is ~ t~'^''^ . The universe shrinks to zero size as t ^ oo. We can also 
calculate expectation value of torsion. In this model the non- vanishing torsion operator is given 
by 

. d 



' dy' 



and one obtains 



(T) 



£(|) 
r(|) 



4\ 3 f3~iAt 



(3.13) 



(3.14) 



.3; (/32+ ^2^2)3 

which shows that {T ) r^ t for small \t\ and |( T )| ~ |tp^ for large |t|. Since the torsion gives 
kinetic energy, the energy density stored in this universe is not zero. In fact, it can be calculated 
as 



<'> = «?■ 



(3.15) 



We have mentioned that the relation to the conventional torsion free gravity is obtained by 
imposing torsion free condition to the wave function, say T^ = 0. However, it is easy to see that 



it is only satisfied by trivial wave functions thus the torsion free condition gives trivial solution 
in this model. 

It is also interesting to see dispersion of an operator A, given by 



A 



A 



{{A-{A)f) = {A')-{A)' 



This quantity measures the quantum fluctuation. We find the dispersion of y™' is 

2" 



a: 



y" 



5 I 2m, ^ 
3 "•" 3 ^ 



r(| + 



r(| + f) 



m 



r(f; 



2+2^ 



and that of the torsion is 



A't 



5 AH' 



r(|) 



V 



l-Hf)Hl 



9{f3^ + AH 



2P 
9(/32 + AH^) 

~ 0.9667. 



(3.16) 



(3.17) 



Then Aym for m < becomes small and Ay grows as |t| -^ cx). Thus in the large \t\ region, 
quantum effect for the y^/tortion is smaller /larger than that in the |t| ~ region. 

We comment that time dependence of a metric similar to (j3.12p have been found in semi- 
classical analysis of the two dimensional gravity coupled to a scalar field proposed by Jackiw 

M, 

S= [ (fx^/^(i){R-K), (3.18) 



where R and A are Ricci scalar and cosmological constant respectably, and (j) is the scalar field. 
Equations of motion become 



i?- A = 0, 

V^^ - A</) = 0. 

Provided ansatz for the metric and the scalar field 

ds"^ = -dt^ + a?{t)dx^, 
(t> = <P{t), 

where a{t) is assumed to take the form o?'{t) ~ t~" with a constant. Then we have 

'(j) - ar^ (j) + Kcj) = 0. 

This equation has the following solution 

(f){t) = t^jj_c+i(VXt), 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



where Ju{z) is the Bessel function. For large t the scalar field increases as (/> ~ t"'^ when a is 
positive. Flat universe corresponds to the case of the A = 0, In this case the scalar field behaves 
(j) r^ i("+i) for large t. In both cases of zero/non-zero cosmological constant, the scale factor a{t) 



10 



shrinks but the scalar field (/)(i) explodes as time evolves when a is positive. Hence the feature 
of this simple model is somehow similar to that of the minisuperspace model discussed in this 
section. The scalar field plays a role of a trigger to the non trivial dynamics in this model. This 
may serve as a classical counterpart for late time behavior of our minisuperspace model, where 
the scalar field is associated to the degree of freedom of the torsion. 

We would like to discuss a possible interpretation of the time dependent wave function 
(|3.10p . Wave function in quantum mechanics is regarded as a one-particle state in the Fock 
space of quantum field theory. According to this conventional interpretation we identify the 
wave function (I3.10p with a one-universe state in the Fock space of matrix quantum gravity. The 
Fock space describes multi-universes, and there are all kinds of one-universe states determined 
by the choices of C(e) in ()3.8p . A time dependent inner product (t'|t) = /q ^(t', y)*^(i, y)dy 
will give a transition probability such that the universe is observed again after a time interval 
t' — t. For the wave function (jS.lOp we easily calculated the transition probability by using (j3.9p 
and found (t'|i) oc (t' — t)~^''^^ . It suggests that this one-universe state is unstable and after 
some time it will decay into another one- universe state with some branching ratio, if we ever 
had introduced operators for creation and annihilation of universes. 

4 Discussion and summary 

In this paper, we proposed a new interpretation of the BFSS matrix model. We applied the HKK 
interpretation, which is originally used for IKKT model by the authors of [9] , to the BFSS matrix 
model. In our interpretation, the time dependent matrix in the BFSS model can be regarded as 
a covariant derivative, and further it is decomposed into geometrical quantities such as vielbein 
or spin connection with explicit time dependence. Therefore, our Hamiltonian determines the 
time evolution of space time. Using this Hamiltonian, we wrote down the Schrodinger equation 
following to the usual recipe for quantum mechanics. 

Several questions for our interpretation of the BFSS model still remain open for further 
investigation. Treatment of the large N limit is one of them. This requires certain regularization 
procedure which is left unspecified in our paper. Such regularization is important not only for 
making more sense out of the equation (|2.2ip . but also for quantitative argument. It is well 
recognized that a derivative operator can only be expressed as a matrix of infinite size. Therefore, 
as was done in this paper, most natural way is to start with infinite A^ theory, and employ a 
regularization scheme later to make matrix size finite. The heat kernel regularization is a possible 
candidate for this purpose. On the contrary, one can also choose an opposite direction — starting 
from finite A^ matrix and taking large A limit later. In this case, we should find a "finite A^ 
version" of covariant derivative which only coincides with true covariant derivative after taking 
large A limit. We notice that the DO-brane field theory developed in y|^ could be useful for 
this problem. 

Relationship between M-theory and our formalism is another issue which is not addressed 
in this paper. Since BFSS matrix model is conjectured to be M-theory, our result of quantum 
gravity would also be understood within M-theory in eleven dimension. We hope to understand 
this issue in future. 
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Another point is a relationship between the Wheeler-DeWitt equation and our Shrodinger 
equation. In subsection 12.31 we have seen that our Shrodinger equation can be regarded as 
a Hamiltonian constraint by introducing a redundant degree of freedom which corresponds to 
rescahng of time coordinate. However, we can stih ask whether our constraint coincides with 
the Wheeler-DeWitt equation derived from a d+1 dimensional theory with manifest covariance 
in whole d + 1 dimensions. Since the BFSS model is obviously non-relativistic, we don't expect 
a direct connection to a covariant theory, at least straightforwardly. Another wishful thinking 
would be to consider a canonical formalism for the d-\- 1 dimensional IKKT model after using the 
HKK interpretation, which hopefully results to a large N effective theory with d+1 dimensional 
covariance. For the moment, we leave this issue as an open question. 

We have also investigated a toy model of universe in detail via a two dimensional minisuper- 
space model. We have observed that our wave function is completely regular and normalizable 
even though the classical metric have a singularity at y ~ oo. Such phenomena of singularity res- 
olution is rather typical in quantum gravity. There exist at least two notions of singularity: the 
singularity appears in the minisuperspace metric and that appears in the Schrodinger equation 
of spacetime (the Wheeler-DeWitt equation) . The singularity which appears in the metric cor- 
responds to either curvature singularity or torsion singularit}o and Einstein's classical relativity 
breaks down, while the singularity which appears in the Schrodiner equation is nonessential but 
simply a regular singular point in an ODE. The regularity can be verified in our wavefunction 
p.lOP , and singularity avoidance is achieved via an exponential fall off of the wave function while 
approaching the point of classical singularity. Similar scenario of singularity avoidance was also 
found in the approach of loop quantum gravity[T6]. Even if we find no regular solution to the 
Schrodiner equation, it is very likely the situation can be improved after including (infinitely) 
many higher spin fields in the large N matrix, which are ignored through this paper for the 
simplification of discussion. 

More comments are in order. Firstly, we can study more realistic setup of higher dimensional 
universe with nonzero cosmological constant. In such case, additional terms which are absent in 
two dimensional case appear in the Hamiltonian. Then it is interesting to make a comparison 
with the wave function of four dimensional universe in ordinary Einstein gravity as discussed 
in [T7] and hopefully we are able to make contact with the cosmological problem by matrix 
model pis]. Secondly, to make a complete and well defined theory in the HKK interpretation, 
we include finite/infinite numbers of higher spin fields eventually [19|. Thirdly, the matrix model 
in the section [3] is nothing but the c = 1 matrix model without potential. One expect that there 
are some hints in this direction as discussed in [20]. Finally, one may extend our study not only 
to a cosmological spacetime, but also to minisuperspace models with black hole geometry. We 
hope to report some results elsewhere. 



To be specific, the curvature singularity happens at y -^ oo if we consider torsionless gravity, however it is 
possible to have vanishing curvature with divergent torsion if torsion is allowed. 
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A Comment on canonical momentum 

In this Appendix we give a formal argument to define the canonical conjugate momentum 
from a Lagrangian which has been given HKK prescription, instead of directly applying HKK 
prescription to the conjugate momentum, To this end it is instructive to work out the procedure 
of quantization in some detail. As we employ the matrix notation (j2.4p . the Lagrangian can be 
written in terms of the adjoint fields as 

1^ — -n 2;(j)^x^ — -/ ^j pti X(^i-)AX[j)BXQ Xjj , l^-J-j 

where h = t^f^t^^ and / is the structure constant of SU{N). We have been employing 

Aq = gauge. The conjugate momentum of x^ is defined as 

7r«^ = ^ = /.^^xg, (A.2) 

where vr'*' = /i vr^' and vr^ is defined through IlJ^*^' = tt]^ i^^. The Legendre transformation 
gives us the BFSS Hamiltonian 

U — ^h ^^ -7r{i)B _i_ '^ fAB fCD uEF^ ^ ^(«)^0) / A Q^ 

n - -tlABT^{i)T^^ 1-' ^■' ^ X(^i)AX(j)BXc Xj^, , (A.Jj 

where Has is an inverse matrix of h . The Poisson bracket of Xij\ and II^-') is written as 

{x«,vrg)}^^ = 5«0)/^^^. (A.4) 

Upon quantization the Poisson bracket becomes the canonical commutator, and vr'*-' acts as 
—id/dxu\ on the wave function ^(x). The Schrodinger equation then becomes 

4^(-) = [-\h''iA^-Q^ + \f''''Ef\h--x,,AX,,,sx':^x^^^ (A.5) 
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The explicit form of the BFSS Lagrangian after the apphcation of the HKK prescription is 
L = - /d'^xv^trJ^ViVi- J[V„V,f| (A.6) 



(r,/ai,)(r,/9L) - ^(i^i/'OfeOl^^r"^™) 1^ • (A-7) 



-i(r,/a^)(r,/9L)-i 

Inspired by the formal resemblance of (|2.10p to (j2.4p . we introduce following notation. Let 
TA represents di and Ojfc, and E^ be e/ and uJi^^ , where the index A runs also / and i,j. Thus 
(j2.10p can be written as 

X(,)(t) = i?(,)^(t)r^, (A.8) 

where we have defined Eu\ = Ru\^Ei . However in order to fullfill the analogy, we have to 
make a crucial assumption that ta and E commute with each other. This is only valid if 
E is not a function of the coordinates x since ta contains the derivative with respect to x. 
Only in this limited situation, the analogy would works. Furthermore, it should be noted that 
the decomposition (jA.SP is simply a formal expression and we do not intend that r can be 
expressed as a matrix though its counterpart, those generators of Lie algebra t'^^, can have a 
matrix representation. With these assumptions we also introduce a metric 



Qab = d x^/hd trgTATB, (A.9) 



and its inverse Q , which correspond to h and Hab respectively. We can rewrite the relevant 
kinetic term in the Langrangian by using these notations 

1 

Then we can compute the canonical momentum 



L = ^QABEf.Ef.. (A.IO) 



vr(,)=a vr(,)B = a j^=E^^, (A.ll) 



'(0 
which gives (|2.16p . The Hamiltonian is 



;GABn^)7:f^) (A.12) 



1 
= -( 
2" 

which gives rise to the kinetic term of ()2.19p . 

As we mentioned above this formal argument works only with the assumptions. Without 

these assumptions it is hard to perform canonical formalism starting from the Lagrangian ()A.7p . 

We regard that one of the possible origin of these subtleties is due to the large N formalism. Thus 

the process performed in section [2] to get the canonical momentum (j2.15p and the Hamiltonian 

(j2.19p by applying the HKK prescription directly to the canonical momentum (j2.5p and the 

Hamiltonian (j2.6p (which is obtained from the Lagrangian of BFSS model (j2.3p ) can be regarded 

as a finite N regularization. 
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B Embedding BFSS into IKKT 

In this Appendix, we demonstrate that the BFSS model can be embedded into the IKKT model 
with the HKK interpretation. It will become clear to us which kind of spacetime dynamics are 
described by the BFSS model with the HKK interpretation. Let us examine only the bosonic 
part of the IKKT action with Lorentzian signature. It is straightforward to include the fermionic 
sector. The IKKT action is 

SiKKT = ^rr[X(o),X(,)][xW,X«] + irr[X(,),X(,)][X«,x(j-)], (B.13) 

where the indices (i), (j) run from 1 to d. 

We consider the following identification in the large N limit: 

X(o) = i<5(o)°Do, Tr= fdt tr, (B.14) 

where Dq = dt + iAq is the gauge covariant derivative and tr is the trace operation of the BFSS 
model. This is nothing but the Lorentzian version of the prescription described by [21] [22]. As 
a result, we have derived the bosonic part of the BFSS action. 

On the other hand, using the HKK interpretation [10], say Xr^\ = iR/^\°''Va, the action 
()B.13P becomes 

S ='^Tr[Va,Vb][V''M] = 1 J dg I d''+'xy^h:^{x,g\iRat"'0,df\x,g). (B.15) 

We have used the relation TrF — > J dg J d'^~^^Xy^Jid+]~{x,g\F\x,g), and ignored torsion here. 
The integral J dg{g\F\g) means the trace operation over the representations of G. We require 
that applying the HKK interpretation to the BFSS action (|B.14p . say Xu-j = iRu\°'^ a^ should 
be consistent with (jB.lSp . This consistency requires the following restrictions : 

i?(0)° = l, i?(o)* = 0, %)° = 0, i?(,/ G 5pm(d), (B.16) 

ei° = l, e/ = 0, e/ = 0, e/e/5i, = /irf(t,x)/j, (B.17) 

ujQ^^ = 0, ujj^^ = 0, (B.18) 

trF = dg d'^x^/h'd {x,g\F\x,g) = i d'^x^/h'dtr^gF. (B.19) 

It is obvious that under these restrictions, we have explicitly broken G = Spin(d, 1) -^ Spin{d) x 
R, where t G R and g e G = Spin{d). We claim that (|B.16P - ()B.19P are the very HKK 
interpretation for the BFSS model. We remark that the second condition in (|B.18p comes from 
the fact that matrix indices (0) and (i) are not mixed with each other in the original BFSS 
Lagrangian therefore it should also be true for the local Lorentz indices and i after using the 
HKK interpretation. The conditions (JB.ITP and (JB.IOP tell us that the HKK interpretation of 
the BFSS model describes a class of curved spacetime equipped with a metric ds^ = —dt'^ + 
/irf(t, x)ijdx^dx^ (/, J = 1, • • • , d). Then the time parameter t in the BFSS quantum mechanics 
indeed corresponds to the time coordinate in the curved spacetime. 
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C Torsion, Curvature 

In this Appendix, we calculate the commutator [Va,Vfe] explicitly and give definitions for the 
curvature tensors and the torsion. 

Recall that the spin connection is introduced in the covariant derivative such that 

V^y'' = d^V'' + uj^\V\ (C.20) 

where /i is a curved space index and a is a local Lorentz index. On the other hand, we have 

V^K = d^Va - UJ^\Vb = d^Va + UJ^J'Vb. (C.21) 

To compute the commutator, we first note VaV^y^ with local Lorentz indices, that is 

(C.22) 
Then the commutator takes following form 

[Va, Vf,]y^ = Tab^iy^y^") + Tlab^'dy'' (C.23) 

with 

Tab"" = daeb^-dbCa^ + UJab^-LOba^, (C.24) 

T^ab^d = (i^(ib {di^^u^d - duUJf.'^d + ^ffe^u^d " ^u^e^ffd) ■ (C.25) 

Here Tab^ gives rise to torsion and TZab'^d is Riemann curvature tensor. The torsion vanishes if 
the vielbein and the spin connection satisfy the torsion free condition. In that case they are no 
longer independent of each other. It is convenient to write (|C.23P in another way: 

[Va,Vb]V' = Tab^{V^V')+nab'dV'' 

= {dae,^ - dbe^t + u^ab^ - u^ba^) {d^V^ + io^ay") 

= Tab'^df.V' + Rab'dV, (C.26) 

where we have defined 

^ab d = da^b'^d - db^ad + '^ae^b^d ~ ^b'^e'^ad + i'^ab ~ ^b^a)^e^d 

= n^.'d + Tab^co.'d- (C.27) 



16 



D Equations of motion 

In this Appendix we show the derivation of matrix equation of motion in HKK interpretation. 
We first recall 

[Va,[Vb,v,]] = [v„(^b,^Ve + 7^f,/^a/)] 

= {VaTbclVe + n.^Tae'V, + llae'^^Odf) 

= [(VaTfc/) + n^'^Tad'' + n,J]Ve + [Tbc^na/^ + VaTZb/^pdf. (D.28) 

The equation of motion in the BFSS model is 

[D',[Do,X^,)]] + [X^^\[X(^^),X^,)]] = 0, (D.29) 

[X«,[X(,),I)o]] = 0. (D.30) 

Second equation gives Gauss low constraint. Let us consider (1D.29P first. After we take Dq = dt 
gauge, applying HKK identification, ()D.29p becomes 

V, + [V^[V„Vi]]=0, (D.31) 

where V means ^Vj = ^e^^dj + ■§^^i'' Ojk- After using (ID.28p . we obtain 

e, 'di + d)/'^o,fc + v'^m + if^n^Ti^ + n^Vj + [^™"7^-,, + v^TZj.kiP'' = o. (D.32) 

We also defined the d-dimensional Ricci tensor as TZ^ = S^^TZkU'' ■ From ()D.32p we obtain 

+ [^.ki + [V^Tfci^' + v'^'nrTiJ + m]iOjki + Tmi'^n^ki + V^jifci) O^^' = 0. (D.33) 

Here we assume that dj and Oki form a part of independent basis of infinite matrix. Thus the 
equation of motion breaks into two independent part: 

e,^ + {V'Tkine/ + 7]'^TkrTiJ + n/ = 0, (D.34) 

Wifcz - e/e/wjfci + ^™"7^™fc; + V%jiM = 0. (D.35) 

For completeness, we rewrite (|D.34p and ()D.35p in terms of R tensor defined in (jC.27p 

e,^ + {V''nine/ + r]'''Tk^"'TlJ-T/u;)J + R/ = 0, (D.36) 
^^kl - h ^ej^ujjki + Tmi'^iRZki - v""''TpJuOjki) + y' {Rjiki - T.j'^uJmki) = 0. (D.37) 
Similarly, (|D.30p becomes 

^e/S/e/ - e/S/e/ + w,^e/ - w,^e/) dj 

+ [d,^t' - ^^t' + ^^ )^t - ^. W + ^J^'' - ^u^'") ^ki = (D-38) 
where diUJ^ ^^ = e/cJ/w^ ^K Finally, we have equations of motion: 

e^die^ -e^die^ + u^teil -d^J'eil = 0, (D.39) 

d,^t' - ^u^t' + ^. % '' - Co, \uo, '' + u^iOo^ - (.iu^f = 0. (D.40) 
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E Proof of the orthogonality ( 13.71) 



We give a proof of the orthogonality (|3.7p . The left hand side of (j3.7p can be written as 



— lini 

2 a^O 



dxxe "■ ""■ J_^2(^x)J_^2(Ve'x) 



'±±K\^-^)'^±± 



2 2 

where we have multiplied the identity e~" ^ \a=Q = 1- Now we apply the formula 

/'°° _ 2 2 1 _p^+<3^ / pq \ 

I dx e"" ^ xJ^ (px) Ji, {qx) = —^e ~^^ ly \ir^) > 



(E.41) 



(E.42) 



that is valid when Rev > — 1 and \Arg a| < f • Here Iy{z) is the modified Bessel function. We 
have 



— lim 

2 a^O 



1 

2^ 



e ia'^ I_|_2 



ee' 
3 \ 2^ 



Using asymptotic formula for luiz) at large z, 



Uz) 



(|E.43P becomes 



lim 



/27rz' 



(y/?-^)^ 



= e 4a^ 



4av TTVee' 
where we have used S{^/e — ve^) = 2^/e6{e — e'). 



5{e-e'), 



(E.43) 



(E.44) 



(E.45) 
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